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S^ . Abstract 



Stationary black holes of massless supergravity theories are described by certain geodesic curves 
on the target space that is obtained after dimensional reduction over time. When the target space 
is a symmetric coset space we make use of the group-theoretical structure to prove that the second 
order geodesic equations are integrable in the sense of Liouville, by explicitly constructing the 
correct amount of Hamiltonians in involution. This implies that the Hamilton-Jacobi formalism 
can be applied, which proves that all such black hole solutions, including non-extremal solutions, 
possess a description in terms of a (fake) superpotential. Furthermore, we improve the existing 
integration method by the construction of a Lax integration algorithm that integrates the second 
order equations in one step instead of the usual two step procedure. We illustrate this technology 
with a specific example 
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1 Introduction 

The study and construction of black hole solutions in supergravity has a long history [1-5]. 
A prominent role in this history is played by solutions that preserve some supersymmetry. 



Such solutions are easier to construct because the preservation of some fraction of super- 
symmetry leads to first order equations instead of the usual second order equations of 
motion. Furthermore, supersymmetric solutions are protected, to some extent, from quan- 
tum corrections and this allowed to compare supergravity states with open string states, 
and thus to verify the celebrated Bekenstein-Hawking law for the black hole entropy from 
direct microstate counting. 

However, non-supersymmetric solutions are perhaps even more interesting. The lack 
of first order BPS equations, that arise from demanding preservation of supersymmetry, 
implies that we need to find other tools to integrate the second order equations of motion 
in this case. This paper is concerned with providing such tools. Our analysis will be 
restricted to spherically symmetric, asymptotically flat black holes. 

The second order differential equations can be derived from two different, but equiv- 
alent, effective actions and the tools that are developed for studying non-supersymmetric 
solutions depend on the choice of the effective action. Let us briefly summarize these 
matters, assuming the simplest case of stationary, spherically symmetric (single-centered) 
black holes. We shall be mainly dealing with static solutions in four dimensions, although 
our analysis in three dimensions will allow us to extend our analysis by introducing a 
Taub-NUT charge. Furthermore we consider four- dimensional solutions, the extension to 
any other dimension being straightforward. 

The first kind of effective action is obtained, for static solutions, by expressing the 
Maxwell field strengths in terms of the quantized magnetic and electric charges m 1 and 
ej via the respective equations of motion (and Bianchi identities). Then all remaining 
degrees of freedom depend on the radial coordinate r, rendering the effective action one- 
dimensional. The action thus obtained describes the motion of a particle subject to an 
external force field described by the effective, non-positive, potential —V((fi,m,e) [6,7] 

C = U 2 + \G rs ^^ s + e 2U V(<p r \m,e), (1) 

where V > is called the effective black hole potential and e 2U is the time-component 
of the static metric. The r -dependence of V arises from the coupling of the scalars to 
the vectors and where we have used <p r to denote the derivative of the fields with respect 
to the radial variable r. In this framework the main tool to describe non-supersymmetric 
extremal solutions is by mimicking the supersymmetric case. This means that one tries to 
find a function W^(U, <f) r ; m,e), named the "fake superpotential" or "fake central charge", 
such that 2e 2U V = \ (duW^) 2 + drW^W^. If this situation is realized the action can be 
rewritten as a sum (and difference) of squares. Putting the separate squares to zero leads 
to first order flow equations 

f = G rs d s W 4 , U = duW A ■ (2) 

This is nothing else but the familiar Hamilton-Jacobi integration method [8], and W4 is 
the Hamilton's characteristic function associated with the autonomous system (1). Such 
road to the understanding of non-susy extremal solutions has been used in e.g. [9-15] and 
references therein. It has later been realized, starting with [16], that this approach could 



be extended to the non-extremal Reissner-Nordstrom black hole, and it has subsequently 
been generalized to other, more involved, non-extremal solutions [17-19] 1 . 

In this paper we use the second form of the effective action, found in [20], which is 
obtained by dimensional reduction of the four-dimensional action over the timelike direction 
to three dimensions. In three dimensions all vectors can be dualised to scalars such that 
one ends up with a supergravity theory, whose bosonic part of the action is given by 

S = J d 3 x v^(|^3 - \9iMWW) ■ ( 3 ) 

The scalars (f> 1 parametrize a target space with metric gij{4>) of indefinite signature [21-24]. 
One can show that gravity decouples from the sigma model such that the non-trivial part 
of the effective action becomes a pure geodesic problem 

C = \ giJ W . (4) 

The target space in three dimensions (with metric g^) contains the target space in four 
dimensions (with metric G rs ) as a subspace. This description of supergravity solutions as 
geodesic curves goes beyond four-dimensional black holes and works for generic supergrav- 
ity solutions that depend effectively on one direction, like e.g. stationary and cosmological 
p-brane solutions, wormholes, instantons, and so on (see [25] and references therein). 

Similar to the first effective action (1), the Hamilton- Jacobi formalism, if it can be 
applied, leads to a (fake) superpotential W^cf) 1 ) such that the second order differential 
equations can be integrated to first order flow equations 

0* = g^djWz . (5) 

As explained in [8], W^U, <f) r ; m,e) is related to Ws(4> 1 ) in a simple manner. We shall 
mainly be concerned with the latter, which we shall simply denote by W. 

Both approaches (l)nd (4) are equivalent, but the second approach (4) is beneficial since 
it makes more symmetries apparent. Moreover, when the target space in three dimensions 
is a symmetric space (as is often the case in extended supergravities 2 ), one can use group 
theoretical algorithms to integrate the second order equations of motion, see e.g. [20,25-47]. 

These algorithms were first established for the case that the target space is Riemannian. 
This case is relevant for the construction of cosmological solutions in supergravity [34, 35, 
39,42]. Later developments have generalized the integration procedures to target spaces of 
indefinite signature [37,40,41,48], that find applications in describing stationary solutions 
such as black holes. Essentially, these algorithms solve the geodesic equations in a two 
step process. In the first step, one relies on the fact that the geodesic equations can be 
rewritten as a Lax pair equation, i.e. as a first order matrix differential equation for the 
tangent velocities. In the second step, one eventually solves the first order system defined 
by the expression for the tangent vector. 



lr The most general form of the non-extremal flow equations were derived in [19]. 

2 The target spaces in supergravity are symmetric for theories with more than 8 supercharges. An 
interesting subset of theories with 8 or less supercharges exhibits symmetric target spaces as well. 



The fact that the geodesic equations on symmetric spaces can be written in a Lax 
pair form, is a strong indication that these equations are Liouville integrable, since this 
is a necessary condition for Liouville integrability [49]. In a recent paper, [48] the Liou- 
ville integrability of the first order equation for the tangent vector was proven, whereas in 
this paper we establish the Liouville integrability of the full second order geodesic equa- 
tions. This we prove by explicitly constructing the correct number of conserved quantities 
(a.k.a Hamiltonians) that mutually Poisson-commute (a.k.a are in involution). These in- 
volutive conserved quantities are appropriate non trivial rational functions of the Noether 
charges. The explicit form of the involutive hamiltonians was discussed in [48] for the 
case of S £(N, M.) algebras (such construction can be shown to work for a variety of other 
symmetric geometries such as the STU model). Its general form will be presented in a 
forthcoming paper [50] . To our knowledge such constructive proof of Liouville integrability 
for this class of models is not present in the existing literature. 

Liouville integrability implies that there always exists a description in terms of a (fake) 
superpotential. This (fake) superpotential can be lifted to the (fake) superpotential in the 
usual four-dimensional formulation that uses a black hole effective potential [8]. Therefore, 
in this paper we prove that all spherically symmetric solutions of symmetric supergravity 
theories where the scalar manifold is a homogeneous coset allow a (fake) superpotential 3 . 
This result is not only important for the understanding of the differential equations that 
govern black hole solutions, but the existence of a fake superpotential, in the extremal 
case, comes with specific implications for the physics of black holes. In particular there are 
implications for the asymptotic stability of the attractor points (for large and small black 
holes) [51]. There are also implications in the context of holography [52,53]. 

As suggested by Liouville integrability, we find a novel way to directly integrate the 
geodesic equations of motion in a one-step procedure to find the expression for the coset 
element. This is a significant improvement of the existing integration algorithms. 

This paper is organised as follows. In section 2 we present the basic mathematical 
notions of symmetric spaces and geodesic equations in Lax pair form. In section 3 the proof 
of Liouville integrability of the second order geodesic equations is presented. Section 4, on 
the other hand deals with the practical integration procedure of the geodesic equations. 
Finally, in section 5 we illustrate and apply this technology to black hole solutions in simple 
supergravity theories. We conclude with a discussion in section 6. 

2 Geodesic curves on symmetric spaces 

2.1 Symmetric spaces 

Consider a real n-dimensional symmetric space G/H*, where G is a semisimple Lie algebra 
and H* is a maximal subgroup. For pseudo-Riemannian manifolds, H* is a non-compact 
real form of some complex semisimple Lie group, whose compact real form will be denoted 



3 This has been conjectured in [18] 



by H . We denote the Lie algebras of G and H* by g, $)* respectively. The orthogonal 4 
complement to Sj* in g is denoted by ^ 

g = f)*®&. (6) 

This decomposition is called the Cartan decomposition and alternatively it can be defined 
by introducing an involutive automorphism 9, the so-called Cartan involution. This Cartan 
involution acts as follows 

0{S?)=tf, 6>(j*) = -j|. (7) 

As 9 preserves the Lie brackets, we obtain 

[S)\S)*]dS)\ [£,£]cio*, [£*,£] C&. (8) 

In case one works in a real matrix representation, the action of 9 can be expressed in terms 
of a diagonal matrix rj (see e.g. [25]): 

VXGg : 9{X) = -r 1 X T r 1 , V = diag(l p , -l q ) . (9) 

Let us also introduce a different decomposition of g that will be useful throughout this 
article. This is the so-called Iwasawa decomposition of g 

g = Sj*®Solv, (10) 

where Solv is a maximal solvable subalgebra of g. The corresponding decomposition at 
the level of group elements does not hold globally on G/H*, since not all elements g of G 
can be expressed in the form g = s ■ h*, where s G exp(Solv) is a solvable group element 
and h* E H*. 

We denote the generators of the solvable algebra by Ta 5 , the generators of Sj* by H a 
and the generators of ^ by K&- The following relations then hold 

K A = \{T A + V Tlr,), (11) 

H a = < Non- vanishing combinations \{Ta — V^aV) f • (12) 

Choosing a parametrization on G/H* amounts to defining a coset representative L(0) e 
G/H*. It is useful to choose this representative in the so-called solvable gauge, meaning 
that it is obtained by exponentiating a generic element of the solvable algebra Solv of (10) 

L(0) = exp(0 J T J ) , (13) 

where we have used the index i to label the 3d scalar fields and consequently also used 
that index for the solvable generators here. Since, as previously mentioned, the Iwasawa 



For symmetric spaces G is semi-simple and we define 'orthogonal' with respect to the Cartan-Killing 
metric of |j. 

5 Or sometimes by Tj. 



decomposition of G with respect to H* holds only locally in G/H*, the solvable coordinates 
0* do not provide a set of global coordinates but only span a patch of the whole manifold. 
However this is not an issue for black hole solutions as we explain in section 5.3. 

For the time being, we restrict ourselves to the solvable patch of G/H*, which is isomet- 
ric to the solvable group manifold exp(Solv), the metric being defined below. The vielbein 
and connection 1-forms, d(f> l Vi A and d0* Wj a resp., are defined via the decomposition 

L _1 dL = d<f) i V i A T A = d<t) i V i A K A + d(f) i W i a H a . (14) 

The metric qab on the tangent space to the manifold at the origin is proportional to the 
restriction of the Cartan-Killing metric on q to & 

9ij = 9abV a Vj B , gAB = aflAB , f] AB = Tr(K A K B ) , (15) 

where V A is the vielbein matrix. The normalization factor a depends on the representation 
of the generators K A and is chosen so as to have the standard normalization for the kinetic 
terms of the scalars in the Lagrangian. The algebraic structure of Solv is encoded in its 
structure constants Jab an d can be either described by the Maurer-Cartan equations for 
the 1-forms V A = dcf) 1 V A or by the commutation relations among the generators T A 

dV A = -y BC A V B AV c , [T A ,T B ] = f AB c T c . (16) 

In view of the following analysis it is useful to define the adjoint representation on Solv of 
the coset representative L(0 J ) e exp(Solv) and of each single generator Ta 

L-'TaL = h A B T B , (T A ) B C = -f AB C ■ (17) 

Then the first equality in (14) can be written in the adjoint representation as follows 

(h- 1 ) B D d l h D c = -V A f AB C , (18) 

where <% = -J^. 

2.2 Geodesic Lax equations 

Let us now discuss the general features of geodesies on G/H* and their description in 
terms of a Lax pair equation. Consider a geodesic on G/H* defined by n-functions (f> l (r) 
of an affine parameter r. In analogy with the Riemannian case [43], it was established 
in [40,41,48] that one can rewrite the geodesic equations on a symmetric space of indefinite 
signature as a matrix differential equation of the Lax form 

±L=[L,W]. (19) 

The Lax operator L and Lax connection W are defined in terms of the pull-back on the 
geodesic of the left-invariant 1-form on G/H*, which we denote by Q 

Q = h- 1 — L = i L" 1 T — h = Y A T A = W + L, (20) 

dr d(j) z 

L = Y A K A , Y A = f ] AB Tr(nK B ), (21) 



where fj AB denotes the inverse of tjab- From the above definition it follows that the 
components Y A of the Lax matrix L are the pull-back on the geodesic of the vielbein 
1- forms 

Y A = ^V A . (22) 

Note that W and L depend on the scalars as well as on their derivatives with respect 
to t. Then the geodesic action can be written as 

S= fdrC= I dT\g AB Y A Y B = /" dr § Tr(LL) , (23) 

where we have used (15). When one works in the solvable gauge, one can choose a repre- 
sentation for which the Lax connection W is given in terms of L as follows 

W = L> - L< , (24) 

where £>(<) denotes the upper-triangular (resp. lower-triangular) part of L. We are 
therefore interested in solving the following Lax equation 

-^L+[L > -L < ,L}=0. (25) 

Algorithmic methods to achieve this have been devised in the mathematical and physical 
literature [43], [37,40,41,48,54-56]. The integration formulae developed, allow one to 
obtain an explicit solution L bo i(t) of (25), obeying the initial condition L so i(t = 0) = Lo- 
in order to extract the solutions for the scalar fields, one should still solve the system of 
differential equations obtained from 

^Tr (Wr))- 1 i- L (0(r))tf A ) f} AB K B = L sol (r) . (26) 

It turns out that, thanks to the use of the solvable gauge, this system of differential 
equations can be solved iteratively [34]. 

The fact that the geodesic equations can be written in a Lax pair form is therefore a 
strong hint of the integrability of the black hole equations of motion, viewed as geodesic 
equations. We should however keep in mind that the Lax equations are first order, while 
the geodesic equations are second order. Only the first integration step has been shown 
to be integrable [48] . The full integration process of the second order equations involves 
two steps and integrability of the first step does not yet establish integrability of the full 
system. 

In section 3, we give a proof of Liouville integrability of the full second order geodesic 
equations. In section 4, we give a novel, more practical, way to obtain solutions of the 
geodesic equations. 



2.3 Noether charges 

Consider the following matrix 

Q = L(r)L(r)L(r)- 1 . (27) 

Using the Lax equation one can infer that Q is a constant of motion: -r* = 0. In fact, 
Q is the matrix of Noether charges and plays an important role in establishing Liouville 
integrability. While L is an element of .ft, Q is an element of the whole algebra g. 

There is an alternative method to describe the geodesic action in terms of the symmetric 
coset matrix 

M(0) = L(0)??L(0) T , (28) 

where r] was introduced in (9). The geodesic action then reads 

S = -- f dr Tr ( — M— M- 1 ) . (29) 

8 J \dr dr J K J 

In this language the solution is immediate (see e.g. [25]) 

M(t) = M(0*(r)) = M(0)e 2QTr , (30) 

where M(0) is the matrix M computed on the values 4> l (0) of the scalar fields 0* at r = 0. 
The geodesic is thus uniquely defined by the matrix Q and the initial point 0*(O) on the 
manifold. From definition (27) it follows that M(0) Q T = QM(0). On the other hand, 
by definition, M(r) should be symmetric for all r and the considered solution implements 
such property. When the initial condition is chosen at the origin, M(0) = i], we find that 

Q T = vQv- (31) 

For generic initial conditions Q is an element of a subspace of q which is isomorphic to ft. 
In what follows we choose the moduli such that M(0) = rj. 

Equation (30) is an explicit integration of the second order equations, however, not a 
useful one for matrices with dimension larger than two as the problem of extracting the 
individual scalar fields from the solution for M becomes non-trivial. Below we present 
an integration procedure that gives the explicit solution for the coset representative L. 
This can then be used to find the profile of the individual scalars. The latter procedure is 
simpler than extracting the scalars out of M because of the upper-triangular structure of L. 
The integration procedure we present below for L is related to the known Lax integration 
algorithms. 

Fortunately, not all questions of black hole physics require the exact profile of the various 
scalar fields such that one can learn already various things in this framework. For instance, 
it is useful to understand the space of black hole solutions and its various subspaces, such 
as 1) the space of all regular solutions, 2) the space of supersymmetric solutions, 3) the 
space of extremal, regular non-supersymmetric solutions, and so on. It turns out that these 
spaces can all be characterized by simple constraints on Q as we explain in section 5. 

9 



3 Integrability 

3.1 Preliminary facts 

Consider the geodesic Lagrangian defined in (23) 

C=\g AB Y A Y B = \g t ^^. (32) 

Note that the above Lagrangian describes an autonomous system in which the radial 
variable r plays the role of time. We define the momenta P, conjugate to <fi l and the 
Hamiltonian function 

dC 
Pi = _ r = ^.^, H = ±g ij P i P j = £. (33) 

d(p l 

The variables 0*, Pj span the phase space of the system and will be collectively denoted 
Z = {0 J , Pj}. The Poisson bracket on phase space is defined as usual 

{0\ <fP] = , {P, Pj} = {P, <fP} = -8J . (34) 

We can write the geodesic equations of motion in a compact form as follows 

Z + {U,Z} = 0. (35) 

From their definition (22) we derive the expression of Y A in terms of the conjugate variables 

Y A = g AB V B i P i , (36) 

where Va denotes the inverse vielbein matrix (Va Vb° g AB = <? y ) of the solvable group 
manifold. From (34) and (36) and the Maurer-Cartan equations (16) one verifies the 
following Poisson brackets 

{Y A , Y b } = -f AB C Yc , { W), Ya} = VA i Qfi, (37) 

where we have defined Ya = gAB Y B , which can be viewed as a basis for the dual solvable 
Lie algebra Solv* . This is the natural Poisson bracket induced by the algebra 6 . 
Let us now define the following n components of the Noether charge matrix Q 

Q A = aTi(QT A ). (38) 



6 For any given Lie algebra g we can turn the dual Lie algebra g* into a Poisson manifold since there 
exists a natural way to define the Poisson bracket. Consider two functions F\, F2 on g* and take Y E g*, 
then we define the bracket {., .} as {Fi, F 2 }{Y) = {d A F 1 ){d B F 2 )!ab c Y c , with d A F = ^-. This Poisson 
bracket indeed reproduces (37) if the representation is taken with the minus sign Ta — > —Ta- This was 
used in [37] to prove Liouvillc integrability of the first order problem for the geodesic tangent velocity 
vector. 
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From the definition (27) of Q one can derive the following relation between Qa and Ya 

Q A = h A B Y B . (39) 

Aside from being constants of motion, the above relation and eq. (37) imply that Qa are 
in involution with the Y A 

{Qa,Y b } = 0, (40) 

and that, moreover, they satisfy the following Poisson relations 

{Qa, Q B } = f A B C Qc. (41) 

Hence, whereas the Ya are related to the Killing vectors associated with the right action 
of Solv, the Qa are related to the left action [57]. 

The Ya and Qa introduced above depend on the phase space variables Z via their 
definitions. In the next section, we show that one can construct a number of constants of 
motion that depend on Ya and Qa- These constants Poisson-commute and one can find 
the correct number of them, so as to constitute a complete set of Hamiltonians as required 
by Liouville integrability. 

3.2 The proof of Liouville integrability 

Liouville integrability, see for instance [58] , is the statement that there exist n functionally 
independent constants of motion 7ii(Z), here referred to as Hamiltonians, that Poisson- 
commute with each other 

{Hi,Hj} = 0. (42) 

This statement implies the usual Hamilton-Jacobi formulation of integrability in terms of 
a (fake) superpotential W, which is the language used in the supergravity literature. This 
is explained below in subsection 3.3. 

What has been established sofar in [48] is the Liouville integrability of the first order 
problem 

Y A + {n,Y A } = 0. (43) 

Since the Poisson bracket on the dual Lie algebra Solv* (the space spanned by the Ya) 
is degenerate, integrability implies the existence of a symplectic foliation for which the 
Hamiltonian flow is Liouville integrable on the symplectic leaves. Each leaf is nothing but 
the co- adjoint orbit of an element (Ya) of Solv*. Denoting the dimension of the coset by 
n and the dimension of the leaves by 2ho (since the symplectic leaves are always even 
dimensional), we have by definition 

2h Q = mnk(f AB C Y c ). (44) 

Reference [48] showed the existence of (n — ho) constants of motion in involution, where ho 
of them correspond to Hamiltonians in involution on the symplectic leaf and the remaining 
n — 2ho constants are referred to as Casimirs. A Casimir %{Ya) is defined by the property 

{U(Y A ),Y B } = , V£? = l,...,n. (45) 

11 



The Casimirs define the foliation since the symplectic leaves are labeled by the values of 
a maximal set of functionally independent Casimir functions. For the moment we are not 
interested in the exact form and construction of these Hamiltonians; we simply assume 
their existence. The explicit expression for the Hamiltonians can be found in [48] and we 
recall the relevant formulas in section 5, when we need them explicitly. 

Let us denote the Hamiltonians on the leaves by H a (Y A ), where a — 1, . . . , ho and the 
Casimirs by Ht{Y A \ where £ = 1, . . . ,n — 2ho- If we use the identities (40) and (41) we 
find 2(n — ho) constants of motion, which Poisson-commute 

n a (Y A ), H e (Y A ), H a (Q A ), H e (Q A ) , (46) 

where the Ti a (Q A ), resp. T-l(.{Q A ) are obtained by replacing Y A by Q A in H a {Y A ), T-L(.{Ya)- 
One can show that each %(,{Q A ) is itself a Casimir, namely it satisfies eq. (45) and 
can thus be expressed as a function of the He(Y A ). To prove this, consider a generic 
representative %(,{Qa) of this set. By construction it satisfies the equations obtained by 
replacing Y A — ¥ Q A in (45) 

m t tr> \* C 



OQ; 



(Q E )fAB C Qc = , V £? = !,..., n. (47) 



Using the relation (39) and the invariance of f A B C under the action of L G exp(Solv), 
we conclude that T-Le(Q A ) = He(jL A B Yb), as a function of Yb, satisfy eq. (45), and thus 
correspond to Casimirs. 

The remaining T-L a (QA) are independent of the 7^(1^) and r H.(.{Y A ). This gives us a 
total of 

(n-h ) + h = n, (48) 

Hamiltonians in involution, thereby proving Liouville integrability of the second order 
problem. The Hamiltonians Hi of eq. (42) are thus explicitly constructed as {Hi} = 
{H a (Y),H i (Y),H a (Q)}. 

3.3 From Liouville to Hamilton— Jacobi and (fake) superpoten- 
tials 

In order to fill the gap between the classical mathematical language of Liouville integrability 
and the language adopted in the current supergravity literature on black hole solutions, it 
is convenient to recall a few basic definitions and concepts concerning the momentum map 
on symplectic manifolds. 

Let M.2n be a real even-dimensional manifold endowed with a closed non-degenerate 
two-form: 

Q = tt a/3 (Z)dZ a A dZ? , 
dn = (49) 
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where Z a denote the 2n coordinates in any given patch. The pair (Aizn , f2) is named a 
symplectic manifold. Consider moreover vector- fields on the manifold M.2n 

X = A-«^, (50) 

with the property that they respect the symplectic structure induced by Q. This means 
that the Lie-derivative of Q along X vanishes: 

= £^n = i x dtt + dixtt. (51) 

In the above equation the symbol ix denotes contraction of the given form along the 
mentioned vector field. Vector fields fulfilling eq.(51) are named symplectic. It follows 
from this definition that for any symplectic vector field X the one-form ix^ is closed and 
hence locally exact on any open neighborhood U C M.2n- In other words on any U we can 
construct a function *}3 X which solves the following equation: 

iy$l\u = #$P\u- (52) 

The map from the tensor product of the tangent space TIA with U into the real numbers: 

< p| r) : TU x U^R (53) 

is named the momentum-map. Clearly the non-trivial topology of the symplectic manifold 
reflects itself into the fact that for the same vector field X the momentum-map on the 
intersection of different patches U and W can be related by a non-trivial transition function 

fuw '■ WflW'i->l Namely we have: 

*Px \urw = fuw^x \ucw- (54) 

Given the existence of the momentum map, one can introduce the Poisson bracket of any 
two functions 7 ^x and *Py associated with two different symplectic vector fields X and 
Y by means of the following: 

WxJy} = ft(X,Y). (55) 

Note that eq.(52) can be interpreted in two ways from left to right or vice- versa. Namely, 
given a symplectic vector field X we can construct its momentum-map representation *}3x 
or, given a function F on the manifold Ai we can look for the symplectic vector field Xj7 
such that: 

*Px F = F (56) 

We can define such a construction the inverse momentum map. 



7 actually these are not true functions rather they are sections of the line-bundle defined by the transi- 
tions functions introduced in eq. (54). 
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A fundamental property of the momentum-map is its equivariance which amounts to 
the following equation: 

V X, Y = symplectic vector fields : Y<£ x = $[y,x] (57) 

Having recalled these concepts let us now return to the case of a canonical system endowed 
with Liouville integrability. According to our previous definitions this means that on the 
symplectic manifold M. of dimension In there exist n independent functions 'Hi(Z) in 
involution, as stated in eq.(42). Let us construct the inverse momentum map of such 
functions, namely n symplectic vector fields Xj such that: 

^ Xl = Hi(Z). (58) 

By construction these vector fields commute with each other and furthermore, in force of 
equivariance, we have: 

XtHj = X# Xj - = ^[x, Xj ] = 0. (59) 

This means that the vector fields Xj are tangential to the n-dimensional level set surface 
T,h defined by the following equations: 

Ui(Z) = hi (60) 

where hi are some set of n real numbers. Since the n independent vector fields Xj provide 
a basis of sections of the tangent bundle to the level surface TE^, it follows that the 
symplectic form Q restricted to E^ vanishes. Indeed we have: 

ft(X,,X,-)k = {1U, Hi) = 0. (61) 

If we use a local canonical patch U C M.2n of the symplectic manifold in which 

n = J2 dp * A d< P = d (J2 Pi d< F) ' ( 62 ) 

i i 

we conclude that on the level surface E^ we have: 

dY^P i d<l> i \\s h = 0. (63) 

Being closed, the restriction to E^ of the one-form: 

P = Pi d(j) 1 , (64) 

is locally exact on the same surface. In other words on each open neighborhood WcS/, 
there exists a local function W^ such that: 

p = dW iu) . (65) 

14 



Let us now consider the definition of the level set surface: 

h i = H i (<f>,P), (66) 

Equations (66) allow us, formally, to rewrite the momenta as functions of the <fi, provided 
det(dhi/dPj) ^ and use the canonical coordinates 0* as independent coordinates on the 
level surface E^. Under these conditions, recalling eq.(33) we arrive at the conclusions that 
(locally) we have: 

0* = ^(0)P,(0 fc , /i fc ), (67) 

which is the more standard Hamilton-Jacobi formulation. The W-function W(0 J , hj) is 
nothing but the Hamilton's characteristic function solution to the Hamilton-Jacobi equa- 
tions [8,58,59] 

We can say that a complete solution W(4> 1 , hj) to the above equation exists in any simply 
connected domain of the phase space in which the Jacobian det(dhi/dPj) ^ is non 
vanishing 8 . 

In problems related to the study of spherically symmetric, asymptotically flat black 
holes, we are not interested in the complete solution W(0 J , hj) to the Hamilton-Jacobi 
equation, since regularity of the four-dimensional solution implies severe restrictions on 
the values of the integrals hi, as we are going to illustrate in section 5.6. 

We shall apply this analysis in section 5.8 to a specific model and prove that in the 
domain spanned by the regular black hole solutions, the Jacobian det(dhi/dPj) ^ is 
non-vanishing. The function W will be explicitly derived for specific solutions. 

We would now like to point out that the integration provided by the Noether charge 
construction 

P 1 (<P,Q) = Qa^- 1 b A V 1 b , (69) 

provides, generically, a non-closed momentum one-form 

dP ^ , (70) 

see [19] for some examples of this. This means that the first order integration provided by 
the Noether charge is different from the first order integration provided by the Hamilton- 
Jacobi formalism and therefore the conclusions reached in [19] about the existence of the 
W-function are incorrect. 

In this way we have shown that Liouville integrability is completely equivalent to the 
existence of the so called (fake) superpotential W. On the other hand, as it was pointed 
out in [48] and as it will be further discussed in a next coming paper [50] , all supergravity 



8 This property is stronger than the so called " flow box theorem", see for instance [60], which implies 
the existence of a complete solution to the Hamilton-Jacobi equation in a neighborhood of any point in 
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which m, m^(o,o). 



theories where the scalar manifold is a generic, not necessarily symmetric coset manifold, 
possess Liouville integrability. Hence for all these cases we proved the existence of the fake 
superpotential. Let us also stress that we can choose one of the Hamiltonians, say Hi, 
to coincide with the quadratic Hamiltonian 1-L(Y A ) = \Y A Y B gAB- Its constant value h\ 
will also be denoted by v 2 in the following and, for black hole solutions, it represents the 
extremality parameter. 

4 Lax integration algorithms 

Proving integrability is one thing; constructing the general integral of the geodesic equa- 
tions is a separate issue. In this section, relying on earlier results [37,40,41,48,54-56], 
we show that this general integral can be constructed. All integration algorithms so far 
developed, focus on giving a solution for the Lax operator L(r). As explained under (25), 
this is in principle sufficient, as it allows to obtain the solutions for the scalar fields after 
solving a second system of differential equations that can be solved explicitly. For practical 
reasons, it is however desirable to circumvent this second integration step, by giving an 
integration formula for the coset representative. In this section, we will propose and prove 
such an integration formula. In section 4.1, we will review how the integration for the Lax 
operator is performed. Relying on that result, the formula for the coset representative will 
be discussed in section 4.2. 



4.1 Integration formula for the Lax operator 

Let us repeat here the results of [40]. We will assume that we work in a matrix represen- 
tation of q consisting of N x iV-matrices. The authors of [40] showed that the solution for 
the matrix elements L pq {r) of the Lax operator can be given in closed form, according to 
a very specific integration formula that depends on an initial value L = L(r = 0) for the 
Lax operator. In fact, several equivalent versions of these formulas were given. 
The integration formula for L(t) that we start from is given by 



L(r) = Q(C) L (Q(C)y 1 , 
where the matrix Q(C) has the following elements 

Ci,i(r) ... C M _i(r) (C^t))u \ 



QiAC) = 



and we have defined 



'QiWQi-iiC) 



Det 



Ci,i(r) 



Ci,i-i(r) 



{CKr))i,j 



®i(C) := Det 



C(r) 
Ci,i(r) 



-2tL, 



o 



Ci,i(r) \ 
C M (r) / 



SVr) 



1. 



(71) 



(72) 

(73) 
(74) 
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The formula (71) can be proven by showing that it obeys the Lax equation. In order 
to do this, we start by using Q(C) to define the following triangular matrices 

*>(C) := Q(C)d(r), (75) 

(^(C))- 1 := Q(C) (cHt))- 1 . (76) 

One can then show that X > (C) (as well as its inverse) is upper-triangular. Similarly, it can 
be seen that X<{C) (as well as its inverse) is lower-triangular. Another useful property of 
these matrices is that their diagonal elements are equal 

(*>(C))« = (*<(C))«. (77) 

The proof of (71) then starts by writing the matrix e~ rL ° in terms of X > (C) and X < (C) 

e~ TLo = (Q(C)y 1 *>(C) , (78) 

e~ TLo = #<(C)Q(C). (79) 

By deriving (78) and (79) with respect to r, one can obtain 

d ,„,„v N _i , , ( d tr>\\ i v /r>\\-l 



C(C) ^ (6(C))" 1 = -L{t) - I ^*>(C) ) (^(C))- 1 , (80) 



Q(C) A (Q(C))- 1 = +L(r) + (^(C))- 1 Qjr^C)) 



51^ 



Using these equations, the triangularity properties of X > (C) and X < (C) and (77), one can 
see that 

Q(C) A (Q(C))- 1 = L>(t) - L<(r) = W(t) . (82) 



By deriving (71) with respect to r, one finds 

d 



dr L ^ + 



Q(C) ^ (Q(C)y 1 , L(r) 



0, (83) 



so that with (82) we find that L(t) indeed obeys the Lax equation. 

The formula (71) gives an integration formula for the Lax operator L(t). As explained 
above, in order to obtain the explicit expressions for the scalar fields, one still needs to 
perform an extra integration. We can actually do better and give an integration formula 
for the coset representative L(r). As the latter is solely expressed in terms of the scalar 
fields and no longer in terms of their derivatives, one can use it to easily extract the scalar 
field solutions. 
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4.2 An integration formula for the coset representative 



We now derive an integration formula for the inverse coset representative. Throughout 
this derivation, we will assume that we are working on the coset space sof - T ' This ^ s n °t 
really a restriction, as one can always find values for p and q and a representation of G, 
such that the coset space G/H* can be embedded in 90 f +g ? 



G_ 



^-» 



SO(p,q) 

SL(p + q) 



SO(p,q) ■ 

As according to this embedding the Lax connection is an element of the algebra of SO(p 
we have 

W(r)r](p,q) = - r]{p , q) W T {t) , 

where 



fj{p,q) :-- 









1, 



p + q = N 



Similarly, the Lax operator L(t) obeys 

L(r)r](p,q) = r](p,q) L t (t) . 

We can now use the fact that L and W can be represented according to the formulae 

L(r) = L>(r) + L<(r) - diag(L(r)) 
W(t) = L^-L^t). 

and substituting eqs. (88-89) into eq. (20) we obtain the following relation 



Lf 



_ x d 



(L(r)) = 2L>(r) -diag(L(r)). 



On the other hand, let us note that by summing (subtracting) eqs. (80) and (from) 
one obtains 



e(c)f q-\c) 



\ [(^(C))- 1 (Jpx^c?) - (4pX>(C)\ (XAC))- 1 

-^ [(*<(c)) _1 (-jfX<(c)) + (^ x >^) (^(^r 1 



4) 

,q), 

85) 
86) 

87) 

88) 
89) 

90) 

81) 

91) 
92) 



Comparing these equations with (82) allows us to express the right hand side of (90) as 



2L>(r) - diag(L(r)) = X>{C) — (*>(C)) 

A simple comparison of equations (93) and (90) shows that 

(L(r)-% = (^(C)L(O)- 1 ).. 

/ Ci,i(r) ... Ci,i_i(r 
=Det 



-i 



93) 
94) 



^(^©i-iCC)' 
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C 



Li-!' 



(C(r)L(0)- 1 ) 1J \ 
r) (C(r)L(0)-%- j 



where we have used that from its definition 

/ C hl (r) ... Cw-i(r) d, 



1 
/ S i (C = 3)«_i(C) 



(*>(C))y = ^ ,„^ ^ D^ 



T 



(95) 



\ C,,i(t) ... C M _!(r) C„(r) / 



Note that we have introduced a constant matrix L(0) G gQ ^ +g ? in the solvable parameter- 
ization of the coset. This matrix can be seen as a matrix of integration constants. The for- 
mula (94) describes an explicit general solution for the coset representative L(r) G so f +g ? • 
This solution is parameterized by the initial data Lq and L(0), which gives precisely the 
correct number of integration constants associated to the geodesic equations. 

5 Applications to black holes 

5.1 Black holes as geodesies 

As explained in the introduction, stationary black holes in d + 1 dimensions can be reduced 
over time to Euclidean solutions in d dimensions. From now on, we assume that d — 3, 
but the extension to any dimension is straightforward, see e.g. [25]. 

Consider a four-dimensional supergravity describing scalar fields <f) r , and ny vector 
fields B 1 , I — 1, . . . , ny. We suppose the scalar manifold to be homogeneous-symmetric of 
the form G4/H4, G4 being the isometry group of the manifold and H4 its maximal compact 
subgroup. The action in four dimensions is given by 9 

S 4 = fU*R 4 - lG rs * d0 r A d0 s - \^ u * G 1 AG J + \vljG 1 A G J ) , (96) 

where G 1 = dB 1 , and G rs , hjj, ujj are symmetric matrices that depend on the scalars 0; 
in particular, G and // are required to be positive definite. The group G4 represents the 
on-shell global symmetry group of the theory, once its non-linear action on the scalars (jf 
is supplemented by a linear (symplectic) electric-magnetic duality transformation on the 
vector field strengths and their magnetic duals. 

Let us now restrict ourselves to stationary solutions having an SO (3) spatial isometry, 
which comprise static black holes though we may allow also for a NUT charge. We then 
"reduce" the action over the time-like direction using the following ansatz 

dsl = -e 2U (dt + A KK ) 2 + e~ 2U dsl, 

B 1 = B J + Z\dt + A KK ) , (97) 



9 We use the more compact form notation here. A p-form A p is written in components as A p = 
■7^1...^ dx Ml A ... A dx^p . The hodge star in d dimensions is defined via, -kdx 1 * 1 A ... A dx^ p = 
(p-d)! £i/ i ■■■•• y d- P tJ ' 1 '" >J ' F dx l/1 A ... A dx Vd - p , where e is the totally anti symmetric tensor (not symbol). This 
implies that ** A p = (— l) p< - d ~ p ) +t A p , where t is the number of timelike dimensions in the <i-dimensional 
space. We furthermore have the relations *A p A B p — -kB p A A p = ^jA A , 1 ... Mp _B M1 - Mp * 1. The exterior 
derivative acts as dA p — -\d f _ tg A^ 1 ,,, tlp dx fl0 A dx^ 1 A ... A dx^". 
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where B 1 and Akk are vectors in d = 3 + and U and Z 1 are scalar fields. Since time is not 
periodic this is not a compactification over time, but the difference between a dimensional 
reduction over a compact or a non-compact direction only enters the KK sector, which is 
assumed truncated here. 

When <1Akk 7^ the 3d vector Akk can not be removed using a coordinate transfor- 
mation and this off-diagonal term in the metric then corresponds to Taub-NUT charge. 

After dimensional reduction we dualize the vector fields B 1 and Akk to scalars Zj, a, 
see Appendix A for the details of this derivation, so to end up with an Euclidean three- 
dimensional theory consisting of a sigma model describing n scalar fields <p l coupled to 
gravity. The fields <p l consist of the four-dimensional scalars <p r , two scalars U, a coming 
from the four- dimensional metric, and 2riy scalar fields Z 1 , Zj originating from the four- 
dimensional vector fields and which will be collected in a 2 ny-symplectic vector Z = 
(Z M ) = {Z 1 , Zj). The target space of the sigma model is a homogeneous-symmetric 
pseudo-Riemannian space of the form G/H*, G being its isometry group and H* a non- 
compact maximal subgroup of G. G is the global symmetry group of the three-dimensional 
theory and contains a subgroup of the form SL(2, M.)e x G4, where SL(2, M)e is the Ehlers 
group acting transitively on the scalars U, a. 

The ansatz for the solutions in three dimensions is such that the equations of motion 
for the scalar fields 0*, decouple from the equations for the metric degrees of freedom. To 
make this clear we use the following ansatz 

ds^ = ^( 3 ) m ndx m dx n = exp[4A(r)]dr 2 + exp[2A(r)]d^ ! ft = ft(j) , (98) 

where (x m ) = (r, 6, tp) are the space coordinates and dQ 2 . is the metric on the unit 2-sphere. 
The radial variable r is chosen so that a/ I ^(3) I 9jg\ ^° e ^-independent. Consistently with our 
symmetry requirements all the scalar fields will depend on r only: ft = 4> % {t). The scalar 
field equations of motion can then be derived from the geodesic action (we refer the reader 
to Appendix A for a derivation of this action) 

£ = bij(<f>W = U 2 + \ G rs <j) r s + \e~ 4U (d + Z T CZ) 2 - \e~ 2U Z T M 4 Z , (99) 

where the negative definite, symmetric, symplectic matrix Ai^ is defined in Appendix A. 
The target space metric is indefinite, the negative- signature directions corresponding to 
dZ M . 

One thus finds that r corresponds to the affine parameter along the geodesic and it is a 
re-parametrization of the usual radial coordinate used for spherically symmetric solutions. 
A geodesic that is parametrized by an affine parameter r moves at constant speed, denoted 
v 2 

\ 9ij {<l>)W = v 2 . (fOO) 

Note that v 2 can be zero, positive or negative, due to the indefinite signature of the target 
space. Note that the effective Lagrangian (99) describes an autonomous system in which the 
radial variable r plays the role of time and eq. (100) represents the Hamiltonian constraint. 
We do not consider timelike geodesies (v 2 < 0) as they can be shown to correspond to over- 
extremal (singular) black hole solutions. Not every geodesic with v 2 > corresponds to a 
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regular black hole solution as we explain in some detail below. For black hole solutions v 
is the extremality parameter and is expressed as the product of the black hole temperature 
T and entropy S [4] 

v = 2ST. (101) 

Hence, extremal solutions, for which v — 0, are described by lightlike geodesies and non- 
extremal solutions by spacelike geodesies. 

The solution for the metric is independent of the details of the sigma model and is 
determined by the geodesic velocity only 10 . It can be found to be 

v 2 >0: exp[2A] = , ! - , (102) 

sinn \vt) 

v 2 = 0: exp[2A] = — . (103) 

T 

The non-trivial part in finding the three-dimensional solutions using the ansatz (98) thus 
consists of solving the geodesic equations for the scalar fields 0\ One of the scalars in three 
dimensions is the Kaluza-Klein dilaton U that comes from the dimensional reduction over 
time. The three-dimensional solutions are uplifted to four dimensions using equations (97). 

5.2 The structure of q 

For homogenous sigma models that appear upon time-like dimensional reduction from 
four dimensions [23,25], the adjoint representation of the isometry group G branches with 
respect to its SL(2,M)e x G4 subgroup, upon dualization of all the vector fields to scalar 
fields as follows 

Adj[G] -* (Adj[SL(2,R) £ ],l)©(l,Adj[G 4 ])©(2,R), (104) 

where R is the symplectic representation in which the electric and magnetic charges of the 
four-dimensional theory transform under the duality action of G4. 
The solvable algebra Solv of (10) has a specific structure [25,42] 

Solv = Solv 2 © Solv 4 © Span(T M ) , (105) 

where Solv 2 = Span(if , Ep ) generates a submanifold S£(2, R)/SO(2), f3 being the high- 
est root of q, while S0IV4 generates the scalar manifold G4/H4 of the four- dimensional 
parent theory. The nilpotent generators Tm, M — 1, . . . , 2ny correspond to the positive 
(restricted) roots 7m defined by the property ^m(Ho) = 1. The space Span(TM) transforms 
in the representation R +i with respect to the subgroup G4 x 0(1, 1) of G, the 0(1, 1) factor 
being generated by H . The algebraic structure of Solv is described as follows 
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[Solv 2 , Solv 4 ] = , [H ,E^] = 2E fi 

[Ho, T M ] = T M , [T M , T N ] = Cmn Ep , (106) 



10 This expression is sensitive to numerical factors that depend on the normalisation of the Einstein- 
Hilbcrt term with respect to the normalisation of the scalar kinetic term. 

21 



where Cmn is the symplectic invariant matrix. 

At this point we have partially fixed normalisations through equation (106). If we 
furthermore fix the normalisation of C such that it has only ±1 as entries then there is one 
scaling freedom left; we can rescale Ep if we accordingly rescale the Tm generators. 

Let us now describe the structure of g and of the subspaces .ft, £)* in the Cartan de- 
composition (6). The matrix r\ defining the decomposition through the involution 9 (7) 
has the following intrinsic expression in terms of the Cartan generator H : 9 = (—1) H °. It 
is easy to verify from (106) that 

7]T M ri = -T M , r]T r r] = T r , r}E Po r] = E Po . (107) 

For notational convenience we define the basis of generators {T4} of Solv as follows: {T4} = 
{T ,T r ,TM,T,} where T = ■^ a , T. = Ep and {T r } is a basis of Solv^. From the above 
properties it follows that the generators Ka of the space .ft, defined in terms of T A by eq. 
(11) have the following form 

{K A } = {K ,K r ,K M ,K.} , K = T , K r = l -{T r + T r T ) , K. = 1 -{T. + T. T ) , 

Km = \{T M -T T M ). (108) 

The non-compact generators K r generate the space $4 in the Cartan decomposition of 
the algebra 04 of G4 with respect to its maximal compact subalgebra S)^. g 4 = i^ 4 © ^ 4 . 
The non-compact generators {Kq, K,} generate the space ^2 of the Cartan decomposition 
of the Ehlers algebra sI(2,M)e with respect to its maximal compact subalgebra u(1)e' 
sI(2,M)e = u(1)e®&2- Finally the compact generators K M span a 2ny-dimensional space 

The metric qab at the origin of the manifold is proportional to the restriction to .ft of the 
Cartan-Killing metric on q, and thus it has 2ny negative signature directions corresponding 
to the compact generators Km- The space A transforms, according to equation (8) in a 
linear representation with respect to the adjoint action of H* . In particular we can consider 
the adjoint action on .ft of the maximal compact subgroup H c of H*. This group has the 
general form H c = U(l)# x iJ 4 , if 4 being the maximal compact subgroup of G4. The 
compact symmetry group H 4 , which acts linearly on the central and matter charges in 
four dimensions, is enhanced in the three-dimensional theory by the \J(1)e subgroup of 
the Ehlers group. With respect to the adjoint action of H c , the subspace .ft Rc spanned by 
{Km}, transforms in a representation R c . This is the representation in which the central 
and matter charges in four dimensions transform with respect to H c . For instance, in the 
maximally supersymmetric theory, G4 = E 7 ( 7 ), H4 = SU(8), G = E 8 (s) and H* = SO* (16). 
The group H c is U(8) = U(1) B x SU(8) and R c = 28 +i + 28_ 1} the subscripts being 
the \J(1)e gradings. In this case we can choose complex basis elements Km, labeled by 
antisymmetric couples [ab], a, b = 1, . . . , 8: (K M ) = (K ab , K ab ). 

The Lax operator L is defined to be an element of .ft. Its components along the genera- 
tors Km are the central and matter charges of the four-dimensional theory, which depend 
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on the quantized charges and the scalar fields. For instance in the maximal theory 

L\a Rc = V2 (Z ab K ab + Z ab K ab ) , (109) 

Z a b being the complex central charges of the four-dimensional M = 8 parent theory. It 
is known that a complex 8x8 matrix can be skew-diagonalized by means of a U(8) 
transformation. The 4 real skew-eigenvalues pk of the central charge matrix Z ab define the 
normal form of the four- dimensional central charges with respect to the action of H c . From 
a three-dimensional point of view, this procedure amounts to rotating L\^ R , by means of 
the adjoint action of H c , into a 4-dimensional abelian subalgebra An = {Kk} of ^r c 

4 

3heH c : h- l L\^ c h = V2Y,PkK k . (110) 

fc=i 

We can generalize the above discussion to a generic theory and define the normal form 
of R c with respect to H c as the minimal subspace {normal subspace) of its basis into 
which any of its elements can be rotated by means of an ^-transformation. Let us de- 
note by p the dimension of this subspace. Thus p is the minimal number of charges into 
which the most general set of central and matter charges, in R c , can be reduced by means 
of an //.--transformation. In the maximal theory we have seen that p — 4. This im- 
plies that the most general black hole solution modulo the action of the three-dimensional 
global symmetry group G (the generating or seed solution) is a p-charge solution. Tak- 
ing ^r c as the basis of the R c -representation, it was shown in [25] that its normal sub- 
space &n = {Ke} is spanned by p-commuting generators Kg, £ = 1, . . . ,p, where p is the 
rank of the symmetric Riemannian coset H* / H c : p = rank(^). In the maximal theory, 

p = rank(^) = rank( U( L ' ) = 4. Similarly we can define a space f)n c of non-compact 
generators from the Cartan decomposition of ft* with respect to the lie algebra S) c of H c : 
$)* = fie © £)r c - Just as ^r c , also i^ Rc defines, with respect to the adjoint action of 
H c , the basis of a representation R c . We can then define a normal subspace $)n — {Jt}, 
£ — 1, . . . ,p, of i^R c with respect to the action of H c . The non-compact generators J a of 
the normal subspace, can be chosen together with the compact counterparts K^ in ^r c 
so that, Ke, Je and Hg = [Kg, Jg\ generate p commuting SL(2, M.) subgroups of G. These 
groups define a subspace of G/H* of the form 

SL(2, R )YC (ln) 



so(i,i)y h* 

which consists of the p-fold product of 0IS2 spaces. This space and its generating isometries 
plays a role in the construction of the seed solutions, relative to the action of G, of regular 
black holes [25]. 

5.3 The solvable parametrization and the Noether charges 

What we would like to stress from the above discussion is that the p commuting compact 
generators Ki in the coset generate a p-torus T p inside G/H*. In the case of a 0IS2 space, for 
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instance, p — 1 and the torus is the non-trivial 1-cycle Si of the corresponding hyperboloid. 
A generic element of T p can not in general be brought to the form of a solvable group 
element of expSolv times an element of H*. This can be easily seen in the case of the 
space dS2- The circle S± at fixed global time interpolates between two solvable patches, 
each describing a copy of Einstein's static universe (see [41] for a detailed discussion of 
this issue). In general each 1-cycle of T p interpolates between different solvable patches 
in G/H*. If one solvable patch is described by a coset representative L(0) e exp(Solv), 
the manifold being homogeneous, we can find a compact element T> of G/H* which maps 
the origin of this patch into a point in any other solvable patch, this being described by 
a coset representative of the form T> ■ L(0) e T> ■ exp(Solv). Thus homogeneity of G/H* 
implies that any two solvable patches are isometric images of each other. As pointed out 
in [29], the subset of points covered by the solvable patches is dense in G/H* and the 
points which are at the boundary of the various patches, and thus are not described by 
the solvable parametrization, are those in which e~ u = 0, where U is the scalar entering 
the metric in the ansatz (97), and that parametrizes the Cartan generator T of Solv, see 
below. Solutions crossing the boundary between two solvable patches are characterized 
by the property that e~ u vanishes at some finite r, and are therefore singular. It was 
shown in [20] for geodesies in (IS2, that only the space-like ones cross the boundary of the 
solvable patch, while the other geodesies unfold within the same solvable patch. Physical 
fields belong to a single solvable patch (physical patch) and, as long as we are interested in 
regular solutions only, we can restrict ourselves to the physical patch only, since in doing 
so, we are not loosing any physics. We can then define the scalar manifold as a solvable 
symmetric pseudo-Riemannian manifold obtained by identifying all solvable patches of 
G/H*: 

M sca i = D\G/H*, (112) 

where D is the discrete subgroup in G/H* whose elements map solvable patches into one 
another. 

Let us now define our solvable parametrization. It is convenient to define the coset 
representative L in the solvable gauge as follows 

L(0 J ) = e- aE ^e^ zMTM h 4 ((l) r )e UH \ (113) 

where L 4 is the coset representative of G4/H4. Using the coset representative (113) one can 
calculate the sigma model metric, the Lax operator L and Lax connection W explicitly. 

The structure constants of Solv, in the basis {Tq,T t ,Tm,T,}, have the following non- 
vanishing components 

c , 



[Ta, T b ] - Jab T c ; 
1 

2 



I 

foM = n &M i /o»* — 1 i frs ', frM = —{T r )M ', ImN* = CmN ■ (H4) 
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JLq — 1 ] JL«o — — 7= Z JL*4 iv j U-*o — — 6 o 5 JL r — JL. 



We find, after some calculations, 

-u 

lL r = — v 2 e (Z T r Li4) ; lL r = — e T t -mnZ Z ; 
h M N = e~ u UM N , V = y/2e- 2U C MN Z N ; L/ = e" 2C/ . (115) 

In (115) we have used the d = 4 coset representative L 4 and the So/t^ generators {T r } 
in two different representations: the Sol t> 4 -adjoint representation (L 4r s , (T r ) s f ) and the 
symplectic representation of the electric-magnetic charges (L 4M - Ar , (T r ) M N ) 

1L 4 i r 1L 4 = 1L 4 r Tg ] 1L 4 ±m L4 = L | 4 m -T !/v ? 

Pr, ^s] — ~(T r )s T t ; [T r , T M ] = —(T r ) M T N ; T rMN = (T r ) M Cpn . (116) 

Recalling the definition of the Lax components V" 4 , L _1 ]L = y A T4, we can now derive 
their explicit form in terms of 0* 



r s ; 



F° = 2 J7 ; y* = v^K 
y^ = V2e- U U N M Z N ; 
y = -e- 2C/ (a + Z T CZ), (117) 

where y r s cic/> r is the vielbein of G4/H4 in the solvable coordinates, G rs = V r l V s u 5^. From 
the Lax operator we compute the sigma model metric 

gafiit = a{Y A Y B fi AB ). (118) 

The normalisation coefficient a can be fixed by comparing with the action obtained from 
dimensional reduction, as given by equation (186) in appendix A 

a = . (119) 

27/00 

We must furthermore have that 

VMM = -V00$NM , V»» = ?700 , Vrs = VOO^rs , (120) 

which can always be achieved by rescalings of the generators. If we insist on the normal- 
isations for the structure constants as given by equation (106) we generically loose the 
freedom to rescale the generators such that r)oo = 1- For this reason we explicitly carry 
factors of fjoo around in our expressions. 

We can now write the explicit form of Qa = Im S Yb 

fj^Qo = 2U-e~ 2U Z T M 4 Z + na = 2M ADM , 

ij^Qr = v / 2(L 4 )/ VfftSsa + 2e~ 2U (Z T T r M 4 Z) + T r MN Z M Z N n , 

VooQm = -V2 (e- 2U M 4MN Z N + C MN Z N n) =V2C MN Q N , 
fj^Q. = -n, (121) 
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where Q = (m 1 , ej) is the vector of quantized charges and, as usual we define M. 4 mn = 
L 4M p L 4 jv P > 0. It is useful to express the Noether charges in terms of the conjugate 
momenta P«: Pu, P r , Pm, P»- We know that Pi = g^ft, which explicitly reads (see 
appendix) 

Pu = 2U , P r = G rs (f) s , Pm = — e " A^4mw2 — ^uCmnZ , P, = -^n, . 

We find, after some algebra 

VooQo = Pu + Z M P M + 2aP. , 

rj^Qr = V2(h 4 ):Vf u P t -2P M T rN M Z N , 

q m = c MN P N -Z M P., 

fj^Q. = -IP.- (122) 

5.4 The (Fake) Super potentials W 4 , W 3 

The (fake) superpotential W4 mentioned in the Introduction, is solution to the Hamilton- 
Jacobi equation associated to the autonomous system (1) describing the radial flow of the 
four-dimensional scalar fields and warp factor: 

l - (duW 4 ) 2 + <9 r W 4 <9 r W 4 = 2e 2U V, (123) 

while W3 = W is solution to eq. (68), associated with the three-dimensional Hamiltonian 
system. The relation between the superpotentials W 4 , W3 = W was found in [8] to be: 

W{(f) 1 ) = W 4 {U,<j) r ;m,e) + Z M C MN Q N = W A + Z I m I -Z I e I . (124) 

5.5 The Hamiltonians 

Let us review, in brief, the algorithm introduced in [48] for constructing the Hamiltonians 
in involution associated with the Lax equation for the S £(N, R)-models 11 (a more general 
construction will be given in [50]). 

In order to construct the Hamiltonians, one makes use of the embedding (84). Let us 
introduce two indices 

\N~ 
a = 0, ••• , — 

a = I,--- ,7V -2a. (125) 

We will denote the Hamiltonians H(Ya) by t) aa . The functions t) aa are then obtained by 
the relation 



det {(L - \l)ij : a + l<i<N,l<j<N-a} 



N-2a v r N 



S a0 X N ' 2a + E baaX"- 2 "-" , « = 0, 



o=l 



2 



(126) 



L As anticipated in the Introduction, this algorithm also works for the STU model. 
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where £ a o is the coefficient of the power A^" 2 ". One can check that the coefficient f)oi 
corresponds to Tr(L) and is therefore zero. The Hamiltonians f)o Q are essentially the 
coefficients of the characteristic polynomial of L and are polynomial in the Y^-variables. 
More specifically, f)o2 is quadratic and corresponds to the Hamiltonian function % defined 
in (33). For a^O, the l) aa are rational functions of Ya- 

5.6 Regularity 

It is important to realise that not every geodesic corresponds to a (physically reasonable) 
black hole solution. We already mentioned that the solutions with v 2 < are over-extreme. 
There is a convenient way to parametrise those Q that correspond to regular black hole 
solutions 12 . This is based on a theorem due to Breitenlohner, Gibbons and Maison [20]. 
This theorem states that all regular non-extremal black hole solutions can be transformed 
to the Schwarzschild solution using ^-transformations. This theorem can be rephrased on 
the level of the matrix Q, as shown in [29]. For the Schwarschild solution, U is the only 
scalar that is turned on. Hence we have L = e UH ° and M = f]e 2UH °. Let us denote the 
Noether charge that belongs to the Schwarzschild solution by Q s , then we have 

Q s = 2vH . (127) 

In the adjoint representation we have the peculiar property that H has a 5 grading struc- 
ture. This means that H is proportional to a diagonal matrix with only —2, — 1, 0, +1, +2 
as entries. If we fix the normalisation of the generators such that Hq contains only the 
—2, —1, 0, +1, +2 as entries then 

Ql = 20v 2 Ql-Uv 4 Q s . (128) 

Since Q transforms in the adjoint of G^\ Q — > gQg^ 1 , we find that the above relation 
is fulfilled for all charge matrices that are connected with Q s and hence for all regular 
non-extremal black hole solutions. For E 8 we have the particular property that the adjoint 
representation is also the fundamental representation. But for all other groups encountered 
in supergravity this is not the case and the particular property of these groups is that the 
fundamental representation has a 3-grading of H ; this implies the simpler equation in the 
fundamental 

Ql = Av 2 Q s . (129) 

In sum, we have that regular black holes obey 

Q 5 = 20v 2 Q 3 -QAv^Q. (130) 

in the adjoint and 

Q 3 = Av 2 Q. (131) 



12 Irregular solutions are not considered physical unless the naked singularity is lightlike such that it 
coincides with the horizon (small black holes). These cases are still of interest since one expects that 
derivative corrections change the solution such as to develop a horizon to cloak the singularity. 
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in the fundamental for all algebras of interest but Eg. From this we find the elegant 
result that extremal regular black holes are described by nilpotent matrices for which the 
expansion e^ T terminates after a finite number of steps leading to polynomial expressions 
for the entries of M and hence to simpler expressions for the scalars. It was conjectured 
in [29] that in the limit v — ¥ another regularity condition exists, which can be written in 
terms of weighted Dynkin diagrams, but we do not investigate this further. 

The regularity condition also facilitates the construction of explicit solutions. If one 
constructs explicit solutions using the formalism of the symmetric coset matrix M, then 
regularity implies 

M = M(0)(i + tQ + ±t 2 Q 2 V (132) 

when v — and 

i+ sMpr lQ+ C os H 2vr)-l Q2 y 

when tJ^O. Similarly in the Lax formalism of section 4, where the computations rely on 
the matrix C(r) = e~ 2rL °, which expands in the same way as the matrix M given above. 
On the level of the Hamiltonians the regularity of the solutions becomes straightforward. 
First we observe that the polynomial Hamiltonians are invariant under G. Hence, from 
Mazur's theorem, we deduce that the polynomial Hamiltonians take the same value as they 
do for the Schwarzschild solution. This means that a necessary condition for regularity is 
that all polynomial Hamiltonians, but the quadratic one (H oc v 2 ), are zero. This condition 
is also sufficient. For that we observe that the polynomial Hamiltonians are the coefficients 
of the characteristic polynomial -Pl(A) 

P L {\) : \ N + H 2 (L)\ N ~ 2 + H 3 (L)\ N - 3 + ... + n N . (134) 

If we use the Cayley-Hamilton theorem, Pl{L) = 0, together with the fact that all "Hj = 
with i > 2, we find 

L N + H 2 L N - 2 = 0. (135) 

We cannot multiply this equation with L 3 ~ N to find the regularity condition (131) since L 
is generically not invertible (and especially not when it corresponds to a regular solution). 
However the above equation does imply that the eigenvalues A of L obey 

\ n + U 2 \ n ~ 2 = 0, (136) 

with solutions A = i^/W^, 0. This however implies the 3-grading structure and hence the 
regularity condition (131) 13 . 



13 Thc regularity condition (131) implies that the solution is on the Schwarzschild orbit, although we 
have not yet proven this in that direction. This is however straightforward since the eigenvalues of L have 
to be 0, \fH2 and —y/T-L^. The multiplicities of the eigenvalues is deduced using TrL = and TrL 2 = 27^2 • 
Indeed has multiplicity N — 2 and the other two eigenvalues each have multiplicity 1, showing that it is 
the Schwarzschild solution. 
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5.7 Normal forms 

The above discussion is an example of how normal forms are useful. By definition a normal 
form is the simplest form of Q G g Oft obtained under the adjoint action of H*. Thus, the 
Schwarzschild solution is the normal form of all regular non-extremal solutions. Another 
useful aspect of normal forms comes from the fact that supersymmetry commutes with 
G and hence the study of the supersymmetry properties of the normal form suffices to 
understand the supersymmetry of the general solution. 

In [25] the normal form of Q for a generic geodesic was derived. This is more general 
then the extended theorem of Mazur [20] since it applies also to extremal solutions and to 
singular solutions. 

The theorem of [25] states that the normal form Qn is given by the following element 
of the algebra g: 

Q N = ® P t=1 [sl(2)eso(l,l)} ©, 9 =1 so(l,l)©iVi/. (137) 

The number p is the rank of the coset H*/H c , with H c the maximal compact subgroup 
of H*. Nil is a nilpotent generator that commutes with the first part of the normal 
form. In [25] the focus was on solutions where Nil is absent and the generating solution is 
described by p axion-dilaton pairs and q dilatons. The reason is that these solutions are 
limits of non-extremal solutions with a (complex) diagonalisable initial condition. It was 
also noted that the known attractor solutions in M = 8 supergravity could be understood 
from the normal form without Nil. Regularity and the absence of naked singularities 
implies that none of the p axion-dilaton pairs can have negative velocity squared. Hence 
extremal solutions are obtained by taking the p axion-dilaton pairs to be lightlike and to 
truncate the q decoupled dilaton scalars. The susy properties of the various solutions for 
the J\f = 8 theory (and the related M = 2 STU model) have been investigated using the 
normal form in [25], see also [29-32] for more recent work. 

5.8 The dilatonic black hole 

Consider the Einstein-Maxwell-dilaton action 

(138) 



s = j ' V\g\{^-\(d<p) 2 -\^F 2 



This is probably the easiest theory that can describe black holes with scalar hair. When 
the dilaton coupling a obeys a 2 = 6 then the action (138) is just the circle reduction of pure 
gravity in d = 5 14 . This model can therefore be embedded in many supergravity theories 
(so does the a 2 = 2 theory). The reduction over time leads to the S£(3, H)/ SO(2, 1)- 
coset. The black hole solutions of this theory have been considered by many authors 
before [19,28,33,61]. Here we redo the analysis in our framework because it nicely illustrates 
our formalism. 



14 In the language of the general sigma model treated in the appendix we have G rs = 1, fiu = e 2 ^ 3< < 
and vu = 0. 
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5.8.1 The S q ( ' 2 jj sigma model 

We define the coset element as in equation (113) and consider the following fundamental 
representation for the generators 



Tn 



Tl 




T. 



To 




T, 



1 



2^ 







(139) 



The notation in the above expressions is as before, where M runs over {1,2} and, since 
the solvable algebra in d = 4 is trivial, we have denoted T r = T^. From these solvable 
generators we can construct the coset generators K& via K& = \(Ta + rfT\r\) with r\ given 
by 

/I 0\ 
77 = ( -1 . (140) 

\0 1/ 



This allows us to compute the coset metric t/ab = T^(KaKb) 



7] = diag(?7oo,^,^ii,r7 22 ,?7..) = diagi 



i i 



i 



2' 2' 2' 



(141) 



This coset metric indeed satisfies the normalisations (120). The coset representative (113) 
is explicitly given by L = e - aT « e v2Z T Me 2v20T^ e 2C/To_ j n com ponents this is 



'."♦# 



L 



V 



o 





-2V2 , 

-2%/2 





x/2. 






a + Z l Z 2 

>/5, 



\ 






From this we can compute the Lax operator 



—iv 2 

2 



V 



2 1 



iv 2 
2 1 

2" 1 



lyO 



"VT 



2" 1 

2 



2^ 



with 



(142) 



(143) 



Y° = 2U. 



^20: 



F^v^e-^+^Z 1 , 

y2 

Y* 



V2e~ u -^Z 2 



-iv. 



a + Z 2 Z X - Z X Z 2 ) . 



(144) 
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5.8.2 The solutions 

In order to investigate regular solutions of this model, we first show how the regularity 
condition (129) can be simplified. As the initial Lax operator Lo = L(t = 0) is an element 
of the S £(3) algebra, it obeys the following identity 

L 3 = ±Tr(Ll)L + ±Tr(Ll)l. (145) 

From (129), one thus sees that the condition for having regular solutions is simply 

Tt(Lg) = 0. (146) 

Of course, this is the S £(3)-version of the statement that all polynomial Hamiltonians, but 
the quadratic one, have to vanish for regular solutions (see the discussion around equation 
(134)). Practically this means that the dilaton charge can be written in terms of the mass 
and the electric-magnetic charges, something that was known before [61]. In the extremal 
case one has 

Tr(L 2 ) = => v = 0, (147) 

we then immediately find that L is nilpotent of degree 3 : Lq = 0. This allows only 2 
arbitrary charges since Tt(Lq) = Tt(Lq) = give 2 conditions for 4 charges. 

In order to find solutions using the integration algorithms, we will use a specific 
parametrization for the initial condition L , inspired from [61], that solves the regular- 
ity condition 

/ 



U 



\ 





a 

2 


— 1 /ct 3 — 4i> 2 a 
2 V a ~P 

=±(a + (3) 

1 //3 3 -4t) 2 /3 

2 V fi-c 


^\ 
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H 


/ a 3 —4v 2 a 
1 a-/3 




/ I3 3 -Av 2 13 
1 P-a 

! ) 



:i48) 



We can choose the gauge in which Ls(0) = 1. This implies that U(0) = 1, which is the 
usual coordinate choice, Z M (0) = a(0) = 0, which can be consistently obtained, when 
n = 0, by shifts of the axions. Only the choice 0(0) = is a restriction. However, if we 
use the dilatation symmetry in 4 dimensions: 

0^0 + c, F Tt ^F Tt ^ c , F H ^ iV"^ c . (149) 

we can generate the solutions with 0(0) ^ 0. 

Let us solve first for the non-extremal solutions. The initial Lax operator is given by 
(148). We present the solution for the initial condition 0(0) = 0. In terms of the two 
functions B(t), C{t) 

a (4v 2 -/3 2 ) + /3 (-4v 2 + a/3) cosh(2t; r) + 2v (a - fi) sinh(2t; r) 



B(t) 
C(t) 



Av 2 (a -(3) 
(-4u 2 + a 2 ) f3 + a (4v 2 -a/3) cosh(2w r) + 2v a (-a + /3) sinh(2v r) 

Av 2 {a - (3) ' 

(150) 

31 



the solution reads 

e 2U = (BC)- 1/2 , e 4, = (B^Cy^ . (151) 

In the extremal limit, defined by v = 0, the corresponding solution is given by equation 
(151), where B and C now read 

B <^ i+ ^ + 2(^r 2 - c w =i -"-2(^- (152) 

We can take the limit where either the electric or the magnetic charge vanishes. Then the 
solution simplifies considerably and the regular horizon collapses to coincide with the black 
hole singularity: the solution becomes a small black hole. One can readily check that the 
nilpotency degree of L becomes 2 instead of 3 in this limit. 

With L 3 (0) = 1 we have Q = L . Through equations (121), this allows us to give a 
physical interpretation to a and (5 

Q = l(a-(3) = M, (153) 

Q, = = -±n, (154) 

Ql — 2\ ~3 = ~ ^Qelectric i (155) 



i cn 3 — Av 2 a 1 

V2 — 2 \ ~o = ~/2^-magnetic i (luOj 

Q<t> = 4(a + f3). (157) 

5.8.3 The normal form 

Since the scalar manifold in three dimensions is S£(3, IR)/ SO(2, 1), the normal form of L 
should be given by 

L e[sl(2)e 80(1,1)] ®so(l,l)& Nil. (158) 

For extremal solutions the decoupled 50 (1, 1) should be zero. One verifies that the s *,' 2 L 
part has nilpotency degree 2 in the fundamental. This implies that we need the extra 
nilpotent generator Nil of degree 3 in order to have a regular solution. Indeed, the solution 
without Nil corresponds to the small black hole solution that has either vanishing magnetic 
or electric charge. Furthermore, if one also truncates the lightlike charges in the S L L part, 
such that one is left with the Nil generator, only then one recovers the double extremal 
solution. This is the solution for which the dilaton is constant everywhere since its initial 
value is at the minimum of the black hole potential and its initial velocity is zero. This 
occurs when a = —(3. 
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5.8.4 Hamiltonian viewpoint 

Following the procedure outlined in section 3, we can construct the Hamiltonians in invo- 
lution needed for integrability, two of which depend on the Y 15 

Hi(y) = -\{ Y y + \{y i y + \{y*y - \{Y*y (159) 

MY) = lY°((Y^-(Yr)-^Y^(Yr + (Y 2 ) 2 ) 

+W3^) 3 " M Y ^ Y * ■ ( 16 °) 

The first of these, Hi(Y A ) corresponds to "H = —^Y^Y A and vanishes for extremal solu- 
tions. The second Hamiltonian, 7i2(Y A ), is cubic and proportional to Tt(Lq). It vanishes 
for regular solutions. 

The proof of integrability requires two extra Hamiltonians, that now depend on the 
charges Qa- In principle they are obtained by making the substitution Y A — > Qa in 
T-L\{Y ), 'H.2{Y A \ This however leads to rather complicated expressions. Since we are 
considering vanishing Taub-NUT charge here, we can adopt a more practical approach. 
We simply choose the two electro-magnetic charges as Hamiltonians 

U 3 (Q A ) = Qi, U 4 (Q a ) = Q2- (161) 

These still Poisson-commute with the Y variables and hence with Tii and % 2 - Moreover, 
for vanishing Taub-NUT charge, they also Poisson-commute among themselves since the 
electric and magnetic charges close a Heisenberg algebra with the Taub-NUT charge. For 
vanishing Taub-NUT charge this Heisenberg algebra degenerates and 1-1% and %& commute. 
On-shell, the Hamiltonians take on constant values "Hj(t) = hi. The interpretations of 
these Hamiltonians is now straightforward. The first Hamiltonian hi corresponds to the 
extremality parameter. The second Hamiltonian keeps track of the regularity condition; if 
hi = (/12 7^ 0) the solution is regular (singular). The third and fourth hamiltonian are 
proportional to the magnetic and electric charge. Hence, regular extremal solutions are 
parameterized by two numbers (Qeiectric, Qmagnetic) and regular non-extremal solutions are 
parameterized by three numbers (Q electric, Q magnetic) and h x . 
Recalling the sigma-model 

C = \ 9ii pp = U 2 + |0 2 - \e~ 2U+ ^{Z 1 ) 2 - l e - 2U -^(Z 2 ) 2 , (162) 

where the generalised coordinates are U, 0, Z 1 , Z 2 , the conjugate momenta are given by 

Pu = 2U, P<t> = 4>, P 1 = - e - 2U+VE ^Z 1 , P 2 = - e - 2U - VE( ^Z 2 . (163) 



15 Note that we are restricting ourselves here to the case where the Taub-NUT charge vanishes. In case 
one allows for a non-zero Taub-NUT charge, the procedure of section 3 leads to an extra Hamiltonian, as 
would be required by integrability in that case. This extra Hamiltonian is a rational function of the Y 
and corresponds to a Casimir, in the terminology of section 3.3. For vanishing Taub-NUT charge, this 
extra Hamiltonian is however not present. 
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The four non- vanishing Hamiltonians in (159) and (161) may then be expressed in terms 
of the coordinates and momenta as follows 

14 _ _1 p2 _ 1 P 2 |_ l 2(7-v / 6</>p2 |_ l 2(7+ v / 6</)p2 
rti — 4 -r t/ 2 r ( j ) -r 2 e r x -|- 2 e .r 2 , 

+ l^U + V-^p up 2 _ _^ e 2U + VEt p ^p2 ^ 

n 3 = &P 1} H 4 = ^P 2 . (164) 

We can write the momenta Pa in terms of the coordinates by solving the 4 equations 

h i = Hi(P A ,{U,(f>,Z 1 ,Z 2 }), (165) 

in terms of the Pa 

P A = f A (U,<t>,Z\Z 2 ), (166) 

when the Jacobian J = det ( -^ ) is non- zero. For the example at hand we find 

7 _ 1 _2[/-a/6</> p2 p i 1 „2U-y/E<pp2p , 1 £ ,VE<t>+2U p2 p , 

J - -§e ^i^ + ^Tge ^fV+ge P 2 ^4> + 

1 p V§<t>+2U p2 p v / 6p2p _|_ 1 p3 hfi7"l 

From this explicit expression one can infer that the regularity condition h 2 = can not 
be satisfied when J = 0. Hence, for physical solutions h 2 = 0, it makes sense to define the 
first order equation (166). 

It then turns out that solving explicitly for the Pa in terms of the coordinates is the 
easiest for the extremal solutions hi — 0. The result is 



:i68) 



Pu = leVS^f +e -7BV 3 ) , 

P^ = &e u (e&hl /S ~ e"75 V 3 ) \Af + e^Wf . (169) 

We can explicitly integrate this to find the 4d fake superpotential 

1 2 2 3/9 

W 4 (C/,0) = 2 6^(0), W(<f)) = -(eV6' t 'hl /3 + e~Vv 4 'hl /3 } . (170) 

As expected the fake superpotential factorizes with the black hole warp factor U [19]. 

In the non-extremal case the expressions for P A in terms of the coordinates are page 
filling unless we take either the magnetic or electric charge equal to zero. When the electric 
charge is zero we find 



P v = -iJ e V6<j>+2U h 2_ 16hu ( 1?1 ) 



P^ = -^^ e V6<f>+2U h 2 4 _ 16hu ( 172 ) 
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where hi < 0. The Ad fake superpotential then becomes 

Wa(U, 0) =^e^+ 2U hi - I6/11 + v /Z ^log(-16/ii + 4 v /Z /^\/e^^ +2t/ /il - 16/11' 

- 4 v /Z /^(£7 + a/|0) • (173) 



This indeed coincides with the expressions found in [19] and in the limit h\ — > it collapses 
to the single charge expression for the extremal superpotential (170) 

/I, 
Vy 4 ([/,0) = e t/ eV2 / l4 . (174) 

6 Discussion 

Let us summarize the main results obtained in this paper. We have established new insights 
in the solvability and integrability of the geodesic equations of motion on symmetric coset 
spaces that appear as sigma models of supergravity theories that are reduced over the 
timelike direction. 

Concerning the solvability we have presented a recursive but closed formula for the coset 
representative describing a generic geodesic solution to equation (94). This extends the 
existing results that only supply a general formula for the Lax operator [34,37,40,41,48,54- 
56]. Finding the coset representative still implied solving first order differential equations 
in that approach. It is not always fully appreciated in the literature that the solvability 
was not yet established for the coset representative L. This is due to the fact that the 
solution for the symmetric combination M = hr]h T can be obtained from exponentiating 
a Lie algebra element. However, in order to extract L in the solvable parametrisation, 
which defines the physical fields [62], from this symmetric combination requires solving an 
involved equation for the compensator. This problem is circumvented by the algorithm we 
propose. 

Concerning the issue of formal integrability we have proven in a constructive way that 
the second order equations of motion, i.e. the full autonomous Hamiltonian system, are 
integrable in the sense of Liouville, once the same property has been established for the first 
order Lax pair problem. This means that we have proven the existence of n constants of 
motion that mutually Poisson-commute, where n is the dimension of the symmetric space, 
starting from those associated with the Lax pair problem. This extends the previous result 
which had proven Liouville integrability of the first order differential equations obeyed 
by the tangent to the geodesic. These results solve an open-standing question about 
the existence of a fake superpotential for black hole solutions. The fake superpotential 
description is the supergravity-way of describing Hamilton-Jacobi integrability and since 
Liouville integrability implies the latter we have proven the (local) existence of a fake 
superpotential for all stationary, spherically symmetric, black hole solutions to symmetric 
supergravity theories. The obstruction to a global existence of W is due to the existence 
of points in which the Jacobian det(dhi/dPj) 7^ vanishes. We have proven in an explicit 
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example that this Jacobian is always non-vanishing in the subspace spanned by regular 
solutions, thus ensuring therein a global existence of W. We expect this to be the case for 
generic symmetric geometries, though a general proof is missing. 

We also presented an explicit construction for the n constants of motion leaning on 
earlier obtained results [48]. The constants are divided in constants that depend solely on 
Y, T-L a (Y), and constants that depend solely on Q, H a {Q)i where Y are the Lax components 
and the Q are the Noether charge components in the solvable directions. The constants that 
depend on Y can be divided in two sets : l-L a {Y) = {Hi(Y),'H a (Y)}, where the He(Y) 
correspond to Casimirs. Among the 1-L a (Y), there are a certain number of polynomial 
constants. These have a simple interpretation. One of them corresponds to the geodesic 
velocity, which is proportional to the extremality parameter ST. The remaining polynomial 
T-L a {Y) then keep track of the regularity of the solution: they have to vanish for a regular 
solution. Among the remaining T-[ a (Y), there is one constant that corresponds to the Taub- 
NUT charge. So far we have been unable to give a physical interpretation to the remaining 
constants in the set T-L a (Y). The constants % a {Q) ar e functions of the electro- magnetic 
charges (e I ,mi) and have no easy interpretation. However, for most cases of interest the 
Taub-NUT charge vanishes and then we have shown that one can replace % a {Q) for the 
pure electro-magnetic charges, making the physical interpretation of all constants, but 
the remaining ones in T-L a (Y), clear. We have explicitly worked this out for the dilatonic 
black holes that arise in Kaluza-Klein theories since this is the easiest possible example to 
demonstrate our findings with. This example has only one non-polynomial / H(Y), which is 
the Taub-NUT charge, so that we were unable to learn about the physical meaning of other 
possible non-polynomial Hamiltonians from this example. We anticipate to investigate the 
Hamiltonians in more involved models, such as the STU model, which should make the 
physical interpretation clear. 

Finally our analysis on Liouville integrability can have a bearing on the quantum de- 
scription of black holes along the lines of [13,27]. In these papers the authors developed 
a quantum description of black hole based on the quantization of the radial evolution of 
the scalar fields and the warp factor. The resulting wave function was suggested to have, 
in the semi-classical limit, the following form: 

*(0*) oc exp(iVW)) =exp(iW 4 (U,(f) r )+i(Z I m I -Z^j)) . (175) 

In this picture the n Hamiltonians in involutions, which we have been dealing with in the 
present paper, will provide a complete set of commuting observables, in terms of which to 
completely characterize the black hole state. 

Understanding the physical meaning of these Hamiltonians is therefore an important 
step in order to understand the physics of black holes in supergravity. 
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A The sigma model in d = 3 + 

In this section we closely follow the seminal paper [20]. 

Consider the action (96). The equations of motion for the vectors B 1 are 

d(tx IJ *dB J -is IJ dB J )=0, (176) 

such that, locally, we can introduce the dual potentials Cj 

flu * dB J - v u dB J = dCi . (177) 

The set of field strengths (dB, dC) are not independent and obey the following duality 
relation 

(dc) = CMi * (dCJ ' 

where 

— 1\ . . __ (fi + vfi~ l v vfi~^ 



CM 4 *( ,^1 , (178) 



, n , ■ Mi=[" ' T -i • (179) 

+i o y y fi v f> J 

Let us now reduce the action over the timelike direction using the ansatz (97). The reduced 
action then becomes 

S 3 = I (\ *R 3 - *dU AdU+ \e w *F KK A F KK - \G rs *d0 r A d<p s 

+ Ifiije' 211 * dZ 1 A dZ J - \e 2U fijj * (G 1 + Z 1 F KK ) A (G J + Z J F KK ) 

- uuiG 1 + Z'Fkk) A dZ J ) , (180) 



The vectors Akk and B 1 can be dualised to scalars x an d Zi by adding them as Lagrange 
multipliers to the action that ensure the Bianchi identities 

S'^Ss + xdFKK + ZjdG 1 . (181) 

Varying the action S' 3 with respect to Fkk and G 1 gives the equations of motion 

dZj = -e 2U * fiuiG 1 + Z'Fkk) - ujjdZ 1 , (182) 

d x = \e AU * F KK + Z I dZ 1 . (183) 
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One can verify that the scalars Zj coincide with the timelike components of the dual 
potentials Cj 

C I = C I + Z I (dt + A KK ), (184) 

where Cj is a vector in d = 3. 

Dualisation of the action S3 is equivalent to eliminating Fkk and G 1 from the action 
S' 3 using the two identities (182, 183). If we furthermore make the field redefinition 

2 X = a + Z'Zj , (185) 

(such that a is a symplectic invariant scalar) we find the action 16 

S 3 = f U * R 3 - *dC/ AdU- \G rs * d<p r A d(f) s + \e~ 2U * dZ T A M±dZ 

- \e- 4U * (da + Z T CdZ) A (da + Z T CdZ)) . (186) 

In the main text we have used the notation 

Z = (Z M ) = (Z 1 , Zj) . (187) 

In a compact notation we write 

S 3 = [\*Rs- \9ij{<t>) * # l A d^ . (188) 



B Conventions for charges and mass 

As usual we define charges via Gauss' law, for which we integrate a closed 3-form over the 
spatial dimensions: 



J3 JS 2 



;i89) 



Q I = ei = l dFj =1 / F T , (190) 

n=l !dF KK = l [ F KK , (191) 

J 3 Js 2 

where 

F 7 = dCj = jujj * G J - is u G J . (192) 

We have denoted magnetic charges by m 1 , electric charges by e/ and the Taub-NUT charge 
by n. There is a subtlety in the definition of the electric and magnetic charge in a spacetime 
for which n^O. In that case the integrals that define the Q's are dependent on the radial 
coordinate and hence not constants 17 . The definition that gives constant charges, when 



16 The sigma model metric corresponds to a Lorentzian version of the metric obtained in [42], where the 
similar problem was considered for time-dependent solutions, i.e. Euclidean sigma models. Once explicit 
solutions for the scalar fields have been obtained, one can uplift these to four dimensions using similar 
uplifting formula as developed in [42]. 

17 The killing one-form e t = dt + Aj<K is not hyper surface orthogonal since it is not closed de t = Fkk 7^ 
0. Hence we cannot define a spatial hypersurface orthogonal to the Killing vector. 
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n^O can, for instance, be derived using Noethers theorem and is given in the main text. 
From the formulas of appendix A we find 

Q N = e~ 2U (CMZ) N , (193) 

n = e - 4C/ (d + Z T CZ), (194) 

where we organised the electric and magnetic charges in the symplectic vector 

Q N = (m I ,e I ). (195) 

The mass of the black hole spacetime is defined by the Komar integral 

M = --L f V a ^dS a p , dS Q/3 = -2n [a r m V^ded(f> , (196) 

where C, a = £* = 1 is the timelike killing vector, n a = n l = e~ u is the timelike normal to 
S 2 ^, r 13 = r T = e u ~ 2A is the spacelike (radial) normal to S 2 ^, and finally, y/a = sin9e 2A ~ 2U 
is the determinant of the induced metric on S 2 . We then find 



M = Tl(0) = U(0). (197) 
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